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ABSTRACT 

The detection of 1.37s pulsations from NuSTAR J095551+6940.8 implies the 
existence of an accreting pulsar, which challenges the conventional understanding 
of ultraluminous X-ray source. This kind of sources are proposed to be massive 
pulsars in this paper. Considering the general relativistic effect, stronger gravity 
of massive pulsars could lead to a larger maximum luminosity, scaled as the Ed¬ 
dington luminosity. The pseudo-Newtonian potential is employed to simulate the 
gravitational held in general relativity, and the Eddington luminosity is calculated 
for self-bound stars (quark star and quark-cluster star) and for the Tolman IV so¬ 
lution. It is found that, for a massive pulsar with radius close to the Schwarszchild 
radius, the Eddington luminosity could be as high as 2 x 10^®ergs“^. It is able 
to account for the X-ray luminosity of NuSTAR J095551-I-6940.8 with reasonable 
beaming factor. It is also suggested that massive pulsar-like compact stars could 
form via this super-Eddington phase of ultraluminous X-ray source. 

Subject headings: accretion — dense matter — pulsars: individual (NuSTAR 
J095551-I-6940.8) — stars: magnetar — stars: neutron 
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1. Introduction 


Ultraluminous X-ray sources (ULXs) have X-ray lumi nosity larger than 10^^ ergs~^, 
which exceeds the Eddington limit of stellar-mass objects (IFeng fc Soria! 1201 ll) . They are 
assumed to be accreting black holes, either intermediate-mass black holes accreting at sub- 
Eddington rate, or stellar-mass black holes with super-Eddington a ccretion. Recently, a n 
ULX in M82 was identihed as a pulsar (NuSTAR J095551-I-6940.8) by Bachetti et ah ( 2014 ). 
with an average period of 1.37s and a 2.5-day sinusoidal modulation. The object has period 
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derivative of —2 x 10“^° s and isotropic X-ray luminosity around 10^‘^ergs“^. Considering 
possible anisotropy in accretion of pulsar, one may think that the actual luminosity would be 
several times of lO^^ergs”^, which is still well above the conventional Eddington luminosity 
~ 2 X 10^®ergs“^ for an accreting pulsar with M = 1.4Mq (Mq is the solar mass). 

The critical luminosity of an accreting compact object is determined by the balance of 
gravity and radiation pressure. Both suppressing radiation pressure or increasing gravity 
would result in a larger critical luminosity. To explain the super-Eddington phenemona, 
on one hand, NuSTAR J095551-I-6940.8 is suggested to be an accreting magne tar, which 
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that is focused in this paper. The gravity in general relativity is effectively stronger than 
the Newtonian one in the limit of strong held, which could lead to a much higher Eddington 
luminosity in principle. This general relativistic effect would be signihcant only for pulsar 
with large mass and small radius so that the real radius is approaching the Schwarzchild 
one. Generally, the radius of self-bound star is smaller than that of gravity-bound star, so 
self-bound stars would have larger Eddington luminosity than gravity-bound star with the 
same mass. Additionally, the critical luminosit y of quark-clu ster star could even be higher 


due to high maximum mass (e.g., see Fig. 6 in iLi et al.ll2015h . 


The modihcation of Eddington luminosity due to general relativistic effect is calculated 
in §2, for two classes of solutions to Einstein’s held equation. One is numerical solution 
based on equation of state (Eo S) for quark st ar and quark-cluster star, and the other is the 
analytical Tolman IV solution (lTolmanlll939l) . The results show that massive pulsars could 
have Eddington luminosity as large as 2 x 10^®ergs“^, which might be able to explain the 
X-ray luminosity of NuSTAR J095551-I-6940.8. 

High luminosity implies high accreting rate, and the central object is more probable 
to have large mass during its evolution. This not only explains the ultraluminous X- 


ray pulsar, but also is consistent with the discovery of 2Mq pulsar flDemorest et al.l 12010 


Antoniadis et ahll2ni3h . The timescale for a 1.4Mq pulsar to reach maximum mass by accre¬ 
tion is also estimated, which has the same scale as that of binary evolution. The spin-down 
rate of NuSTAR J0955 51-I-6940.8 m ight be similar to an accreting pulsar with normal mag¬ 
netic field of ~ 10^^ G ( Tong 2015 1. Thus NuSTAR J095551-I-6940.8 could be an accreting 
massive quark star or quark-cluster star, without the need for strong magnetic field in the 
magnetar models. If future observations could find evidence for whether strong magnetic 
field exists or not, it will be possible to distinguish accreting self-bound star and accreting 
magnetar. 


The modihed Eddington luminosity in strong gravity is calculated in §2. In §3, the spin 
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behavior and the timescale of accretion are discnssed. Conclnsions are made in S4. 


2. Eddington luminosity with pseudo-Newtonian potential 
2.1. Modified Eddington luminosity 

A psendo-Newtonian potential 'h = —GM/{r — Rg) (the Schwarzschild radins Rg = 
2GM/c^) is nsed to simnlate th e grav itational field of the Schwarzschild metric, as first 
proposed bv iPaczvhski fc Wiita fjl980f) . The potential together with the Newtonian eqna- 
tion of motion conld have resnlt which is in good agreement with that of general rela¬ 
tivity (ILuI I1985[) . and largely simplify the calcnlation. So it has been widely nsed when 
dealing with black hole accretion disks. The Eddington Inminosity is derived by balanc¬ 
ing the gravity /g and radiation pressnre f\. In the case of fnlly ionized hydrogen plasma, 

/g ~ 


-m, 


,VT = —GMm-p/{r — Rg)^, and /r = pa, where m-p is the proton mass, p is the 
energy density of radiation field, and a is the Thomson cross section ctt with non-relativistic 
electrons and weak magnetic field. From fg + fr = 0, we have p = GMmp/{r — Rg^aT, 
and the maximnm Inminosity (i.e., modified Eddington Inminosity) emitted from the star’s 
snrface, for the case of isostropic accretion, is then 


^, n AirmncGM 
^Edd,s = pc X dvrr = X 

(Tt 


R. 


— LEdd X (1 — Rs/r) 


-2 


( 1 ) 


As a resnlt of gravitational redshift, however, the Inminosity observed at inhnity is. 


-^Edd ~ -^Edd, s X (1 Rs/t) — i^Edd X (1 Rs/t) 


-1 


( 2 ) 


Compared with the conventional Eddington Inminosity LEdd = dvrmpcGM/crT, there is 
a correction factor (1 — Rg/r)~^ in the modihed Eddington Inminosity L^^^. The correction 
factor (1 — Rg/r)~^ increases rap idly as r/Rg app roaches 1. When r/Rg < 9/8, the central 
density of star wonld be infinite (IBnchdahll 119591) . so an npper limit is set to the correction 
factor, (1 — Rg/r)~^ < 9. The correction factor is nnimportant at large r; it has a valne of 2 
when r/Rg = 2, and 1.25 when r/Rg = 5. It is then worth noting that, if a star is extremely 
compact (r is small), its L^^^ could be one order of magnitude larger than LEdd- 


Moreoyer, maximum luminosity approaching or even higher than L^^^ can only be 
possible in case of anisotropic accretion. In case of isotropic accretion, the general luminosity 
cannot be larger than LEdd with which the gravity and radiation can certainly balance at 
sufficiently large distance where Newtonian gravity works (note: LEdd is independent of 
distance r, while L'^^j decreases as r increases). Nevertheless, if the spherical symmetry is 
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broken, radiative photons conld be opaqne in the direction of accretion flow bnt transparent 
in other ways so that the balance of radiation pressnre and gravity wonld be achieved aronnd 
the star snrface where the radiation pressnre would be stronger than gravity at larger radius. 
Inside the accreted fluid, the radiation intensity should decrease with increasing distance. 
The radiation pressure could be weakened by the above factors, and balance of pressure and 
gravity at larger radius is then possible. Therefore, the value of at the radius R of the 
central compact star is chosen to scale the star’s maximum accretion rate. 


The modihcation of Eddington luminosity due to general relativistic effect is discussed 
for mass-radius relations in two cases. One is numerical result based on EoS for quark star 
and quar k-cluster star, and the other is the analytical Tolman IV solution to Einstein’s held 
equation flTolmanlll939l) . 


2.2. Case I: quark-cluster star and quark star 


The EoS of dense matter at supra-nuclear density is still not well known because this key 
question is essenti ally of non-perturbative qu antum chromo-dynamics in the regime of low- 
energy scale fe.g.. ILattimer &: Prakashll2004l) l. Among different models for the i nner struc¬ 
ture of p ulsar, a quark-cluster star could n aturally have very stiff equation of state flLai &: Xu 


2009a b; Lai et ah 20131: Guo et al 
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massive pulsars f Pemorest et al. 201ol: 


which is consistent with the later discoveries of 


Antoniadis et al.ll2013l) . Recently, an analysis of the 


type I X-ray burst observations may show that the compact star in 4U 1746-37 could have 
an extremel y low mass ari d small radius, that could also be understood in the quark-cluster 
star model (iLi et al.ll2015[) . 


Various models for pulsar can be divided into two classes, self-bound and gravity-bound 
stars. Gravity-bound star includes conventional neutron star and hybrid star, th e radius of 


which usually decr e ases as the mass incr eases. Self-bound star in cludes qu ark star fllothlll970 


Alcock et al.lll986l: iHaensel et al.lll986fl and quark-cluster star (jXull2003l) . where the radius 
generally increases with mass. The general relativistic effect would be more important for 
quark star and quark-cluster star, since they usually have smaller radius than gravity-bound 
star of the same mass. The modihed Eddington luminosity is calculated for several EoS of 
self-bound star. 

The hydrostatic equi librium equation in general reP tivity is the Tolman-Oppenheimer- 
Volkoff (TOV) equation fjOppenheimer fc Volkoff 1939 1, 


dP 

dr 


Gm{r)p (1 + ^)(1 + 


1 - 


2Gm{r) 


( 3 ) 
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where m(r) = p-Aur''^dr'. Combine the TOV equation with the EoS, the mass and radius 
of the star can be calculated by numerical integration. For quark star, the simple MIT 
bag model is chosen. For quark-c luster star, Lenn ard-Jones potential flLai &: Xull2009al) and 
a corresponding-state approach flGuo et ahl 120141 1 have been used to derive the EoS. The 
modified Eddington luminosity is calculated for the mass-radius relations given by the 
above EoS, and the results are plotted against mass in Figure [TJ When the mass is small, 


Lgdd is similar for different EoS. The MIT bag constant is chosen to be 57MeVfm“'^ so 
that we have has a maximum mass of ~ ‘2Mq, and the corresponding critical luminosity is 
~ 5 X 10^®ergs“^. Quark-cluster stars have maximum mass in the range of 3 — 4:Mq, and 
Lgdd could be as high as 2 x lO^^ergs”^. 


2.3. Case II: Tolman IV solution 


Besides the numerical solutions employing physically motivated EoS, analytical solutions 
to Einstein’s held equation can also be obtained with appropriate conditions on the metric. 
This mathematically motivated method was hrst suggested by Tolman, and a number of new 
solut ions were obta ined, among which Tolman IV solution is suitable to describe self-bound 
star (jTolmarjll939h . Based on the mass and radius given by the Tolman IV solution, T^dd 
is also calculated, for more general results independent of detailed EoS. 


For static and spherical distribution of matter, the line element can be written in the 
simple form, 

= —e'^dr^ — r^dd^ — sin^ dd0^ -|- e^df^. (4) 

Tolman IV solution stems from the assumption. 


e" v' j^r = const. 


( 5 ) 


and the result is determined by two independent parameters A and C. 
Pc central pressure pc, radius d?, surface density pb and mass M are 

The central density 

3 3, 11 

^ + and SttPc = ^ 

( 6 ) 

= 3(72(1 - 

( 7 ) 

Pb = Pc-5pc + 8p^/(pc + Pc), 

(8) 

R\ {1-RVC?)(1 + RVA^)\ 

2 \ l + 2i?7A2 J’ 

( 9 ) 


respectively (in this subsection, the G = c = 1 units are used). From equation (I7l) and 













6 




Fig. 1.— Upper panel: the mass-radius curves for several EoS. Lower panel: the relation 
between the Eddington lumi nosity and mass M for these EoS. Solid l ines represent Eo S 
from GLX (IGuo et ahl 120141 ). and dashed lines represent EoS from LX (ILai fc Xijl2nn9al ). 
Different numbers indicate different parameters. The dot-dashed line is the result of MIT 
bag model (bag constant B = 57MeV fm“^). 
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equation ([9]), 


2MIR 


{l-R^/C^){l + R^/A^) 

1 + 2i?V7l2 

H(1 - A^C^) < I 


( 10 ) 


So there is a lower limit R/2M > 1.5, and the correction factor (l—Rg/R) ^ = {1 — 2M/R) ^ 
could not be larger than 3. 


Since there are two independent parameters, the value of surface density is hxed to 
reduce one degree of freedom, and thus all the properties can be regarded as functions 
of the single parameter A. For self-bound star with hnite surface density, three cases pb = 
1.5ps, 2ps, 2.5ps are considered, where ps is the nuclear saturation density. The resulting pairs 
of M and R are plotted in Figure [2|, and the modihed Eddington luminosity for different 
mass is also shown. For Tolman IV solution, is smaller than that in the case of quark- 
cluster star. If considering the beaming effect, it might as well explain the X-ray luminosity 
of NuSTAR J095551+6940.8. 


Discussions 


NuSTAR J095551-I-6940.8 has a spin period of P = 1.37s and a period derivative of 
P ^ —2 X10“^°. The pulsar is spinning up and gaining angular momentum from the accretion 
flow. According to the conservation of angular momentum, the strength of dipole magnetic 
hied Bp can be estimated. For a pulsar with mass of 2 Mq and beaming factor of 0.2, the 
value of Bp is 1.6 x G , where Rq = P/10® cm and L 40 is the X-ray luminosity in 

units of 10^®ergs“^ ( Tong 2015 ). So the spin evolution of NuSTAR J095551-I-6940.8 might 
be similar to an accreting pulsar with normal magnetic held of ~ 10^^ G, and there might 
be no need for strong magnetic held being necessary in the magnetar models. 

The birth rate of accreting quark star or quark-cluster star and accreting magnetar would 
be diherent, which might be studied in the future work. This could answer the question, 
what proportion does pulsar system take up in the ULX population? What is the nature of 
ULX pulsar, accreting self-bound star or accreting magnetar? The problem is also related 
to the nature of anomalous X-ray pulsar/soft Gamma- ray repeater: are they magnetar or 
“quark star/quark-cluster star -|- fallback disk system” flTong &: Xull201ll) ? 
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Fig. 2.— Upper panel: the mass M and radins R given by different valnes of parameter A in 
Tolman IV solntion. Lower panel: the relation between the modihed Eddington Inminosity 
and mass for Tolman IV solntion. Different lines represent three cases of snrface density pb 
— l.Sps (dashed line), 2ps (solid line) and 2.5ps (dot-dashed line). 
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3.1. Accretion process and formation of massive pulsars 


From calculations in §2, only massive pulsars above 2.5Mq could have critical luminos¬ 
ity larger than ~ 10^®ergs“^. From the perspective of evolution, the pulsar is also more 
probable to have large mass since high luminosity indicates high accretion rate. When 
born in supernova explosion, the pulsar may probably have a mass of IAMq only. Accretion 
from companion star may lead to the formati on of mass i ve pul sars with mass of 2Mq or even 
higher. According to mass-radius relation in I Guo et ahl (12014J) . the time scale is estimated as 
following for a IAMq pulsar to reach the maximum mass with accretion rate of MEdd, where 
MEdd is the critical accretion rate corresponding to g. The accretion of pulsar might 
be anisotropic because of strong dipole magnetic held (~ 10^^ G), and th e critical luminos¬ 
ity w ould then be larger than the isotropic case due to geometric effect flBasko fc Sunvaev 
19761 ). Gertainly, it could be possible that accre tion toward an compact obje ct occurs at 
super-Eddington luminosity in a binary system fje.g.. Wiktorowicz et al.N2015l) . The situ¬ 
ation of lOiifEdd accretion is also considered, and the results are shown in Figure [31 The 
timescale is ~ 10^ years for MEdd and ~ 10^ years for lOMEddi w hich is consistent with that 
of binary evolution models fjShao fc Lill2015l: iFragos et al.ll2015l) . 


Gonsider a binary system composed of a high-mass star A (> 12Mq) and a low or 
intermediate-mass star B (~ 4 M 0 ). The high-mass star A hrst evolves and explodes to 
become a pulsar with mass of 1.4Mq, and the mass of companion B could be increased 
to ~ 5Mq through accretion When star B expands and exceeds its Roche lobe, mass will 
be transferred to pulsar A via the Lagrangian point. As the mass of pulsar approaches 
2Mq, the accretion rate would become larger and it behaves as an ULX pulsar. Then 
in the common envelope stage, the orbit decays with the release of orbital gravitational 
energy and the envelope could be blown off. Finally there leave a massive pulsar and a 
white dw arf of ^ O.SMrr^. The above scenario mig ht explain both t he discoveries of 2M fn 


pulsar flDemorest et al.ll2010l:lAntoniadis et al.ll2013l) and ULX pulsar (IBachetti et al.ll2014f) . 


4. Conclusion 

The identihcation of NuSTAR J095551-I-6940.8 as a pulsar challenges the understanding 
of accreting neutron star. Gonsidering the general relativistic effect, the Eddington luminos¬ 
ity of massive pulsars could be much larger than 10^® ergs“^ in the frame of self-bound star. 
For numerical solutions based on detailed EoS, the Eddington luminosity of quark-cluster 
star could be as high as 2 x 10^® ergs“^, and that of quark star can reach 5 x 10^® ergs“^. For 
analytical solutions to Einstein’s held equations, the Tolman IV solution gives Eddington 
luminosity up to 10^® erg s“^. The X-ray luminosity of NuSTAR J095551-I-6940.8 can be well 
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Fig. 3.— The growth of mass with time at accretion rate of MEdd and lOMEdd, from 1. 4 
Mq to the maximum mass. Different lines are the results of different EoS flGuo et al.ll2014j) . 
The horizontal line detonates the position of 2Mq. 
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explained by massive pulsar with critical accretion rate. 

In addition, the spin behavior of NuSTAR J095551+6940.8 can be accounted for by 
normal magnetic hied. The timescale for a 1.4M0 pulsar to reach maximum mass by accre¬ 
tion is consistent with that of binary evolution. It is also proposed that as a binary system 
evolves, it might undergo the ULX phase in the accretion process, and dually a massive 
pulsar (> 2Mq) could be formed. 


We would like to thank the pulsar group of PKU for useful discussions. This work 
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